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We present a new and complete analysis of the n-bounce resonance and chaotic scattering in 
solitary wave collisions. In these phenomena, the speed at which a wave exits a collision depends in 
a complicated fractal way on its input speed. We present a new asymptotic analysis of collective- 
coordinate ODEs, reduced models that reproduce the dynamics of these systems. We reduce the 
ODEs to discrete-time iterated separatrix maps and obtain new quantitative results unraveling the 
fractal structure of the scattering behavior. These phenomena have been observed repeatedly in 
many solitary-wave systems over 25 years. 
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Solitary waves — localized disturbances that travel with 
unchanging shape and velocity — are ubiquitous in phys- 
ical science, and are seen, for example, in fluid mechan- 
ics, optics, solid-state electronics, and even quantum field 
theory. A natural question is what happens when the 
wave hits an obstacle or two such waves collide. 

In dissipative systems such as electrical signal propaga- 
tion in nerve fibers or reaction-diffusion systems, two in- 
teracting waves generally merge into a single larger wave. 
In completely integrable, or soliton, equations, by con- 
trast, interacting solitary waves emerge from a collision 
intact and with their original speeds, but a slight shift 
in their position, which is well-understood through the 
theory of inverse scattering. 

Collisions in dispersive wave systems that are nei- 
ther dissipative nor completely integrable may produce a 
much wider range of behaviors. We focus on one, the 2- 
bounce, or, more generally n-bounce phenomenon. Two 
counterpropagating waves with sufficient relative initial 
speed (or one wave incident on a localized defect) will 
pass by or reflect off each other with little interaction, 
while for most initial speeds below some critical velocity 
Vc they will become trapped, forming a localized bound 
state. At certain velocities below Vc, the waves become 
trapped, begin to move apart, and come together a sec- 
ond time before finally moving apart for good — the so- 
called 2-bounce solutions. In addition to the 2-bounce 
resonant solutions, one often finds 3-, 4-, or, more gen- 
erally, n-bounce solutions. Figure shows a 2-bounce 
resonant solution to ([T]), and figure [ip shows the sensitive 
dependence of the final speed on the initial speed, with 
the number of 'bounces' indicated by color. The initial 
conditions leading to these behaviors are interleaved in 
a manner often described as fractal. This was first seen 
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FIG. 1: Color onhne: (a) A '2-bounce' solution to PDE ([T}. 
(b) Vin vs. Vout for kink-antikink collisions to ([T| showing 
chaotic scattering, (c) Two-bounce solution to ODE model 
for ^ {X{t) solid A{t) dashed.) (d) uin vs. Wout for the ODE 
model, annotated as in (b). 



in kink-antikink collisions in the mid 1970's (see [T] and 
references therein), subsequently found in models from 
astrophysics , optical fiber communications [3] and per- 
haps most recently in 2007 in collisions between topolog- 
ical solitons arising in quantum field theory |4|. Figure 
shows a 2-bounce resonant solution of the model ODEs 
for system ([T| (discussed below), and figure [iji shows 
that the ODE model reproduces the fractal interaction 
structure of the PDE, if not the exact structure. 

We analyze these phenomena through systematic 
asymptotics applied to 'collective coordinate' models, 
low-dimensional model systems of ordinary differential 
equations (ODEs) derived from a variational princi- 
ple that reproduce the dynamics in numerical simula- 
tions. We construct, using Melnikov-integrals and for- 
mal matching procedures, approximate n-bounce reso- 
nant solutions to the ODEs and derive an iterated map 
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that explains the fractal structure [5] [B]. 

Previous studies treat the results of numerical simu- 
lations (both ODE and PDE) as experimental data and 
have remarkable success analyzing these results using a 
combination of numerical simulation and ad-hoc calcula- 
tions. They derive approximate resonant velocities using 
least-squares fitting of numerical data. By contrast, we 
obtain formulae dependent only on the equations' param- 
eters, and not on any empirical constants. 

The phenomenon was first observed in kink-antikink 
collisions in nonlinear wave equations by Campbell, 
Peyrard et al. ^j, including the 0^ equation, 
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Figure [TJd is a new computation that reproduces one from 
their first paper. What was in the early 1980's a very dif- 
ficult and time-consuming computation we reproduced in 
a short time on a PC, with improved detail showing nar- 
rower n-bounce windows between the primary 2-bounce 
windows [6.. Fei, Kivshar, and Vazquez subsequently ob- 
served 2-bounce solutions in collisions of kinks with Dirac 
delta potentials in the sine-Gordon and 0"* equations [7] , 

Utt — Uxx + (1 ^ eS{x)) sinu = and (2) 
Utt - + - £S{x)){u - u^) = 0. (3) 

Tan and Yang saw it in collisions between orthogonally- 
polarized solitons in birefringent optical fibers |3j, de- 
scribed by coupled nonlinear Schrodinger equations: 

idtu, + dlu^ + (|u,p + e\ui_i\^)u^ = 0, i = 0, 1. (4) 

What all these (non-integrable) dispersive wave equa- 
tions have in common is a second mode which can draw 
energy from the propagating wave. When the solitary 
wave is taken to model a pseudo-particle, this corre- 
sponds to an internal oscillatory degree-of freedom. This 
transfer creates an effective energy barrier, preventing 
slow waves from escaping the collision location. 

In the following paragraphs, we analyze the behavior 
shown in figures [TJ; and [TJi. The general form of the 
ODE model is given below in ^ . We provide the critical 
velocity for capture in (111 and displayed in figure [sJl 
and[3]3. The locations of the 2-bounce windows, and the 
narrower 3-bounce windows is given in equation ( |19| as 
special solutions of an iterated map we define below. 

Each system above is well-known to possess a varia- 
tional form [S]: their solutions minimize a Lagrangian 



For systems ([T])-(|4|, the ansatz depends on a variable 
X{t), parameterizing the distance between the two inter- 
acting solitary pulses in systems ([T]) and ^ and the pulse 
position in (2]) and (|3|, and a variable ^(i) measuring 
the amplitude of a second mode of oscillation — two such 
modes in the ODE model equations for system ([s]). The 
ODE models take the general form (after some rescalings) 

mX + U'{X)+F'{X)A^Q] A + uj"^ A + cF{X) = (6) 

where c or cj^^ is a small parameter, allowing the use of 
perturbation methods. The ODE model for systems ([ij 
and ([3| contain additional terms but can be treated using 
the methods described herein. We refine our terms de- 
scribing the ODE model: an n-bounce solution is one in 
which X{t) escapes to infinity after n interactions, and 
n-bounce resonance is such a solution for which, addi- 
tionally A(t) — > as i — > ±c» so that = Wout- 

We consider two such models here. After rescaling 
time, system (|2| is modeled by the equations 

U{X) = -2 sech^ X; F{X) = -2 tanh X sech X; (7) 

where m — A, uj"^ — 2/e — e/2 and c — e. The ODEs 
that model equation ([T]) are algebraically complex and 
are studied in [2, 6 , but their essential dynamics (deter- 
mined by the topology of the phase-space) are captured 
by making U {X) the Morse potential and choosing a sim- 
ple F{X) which vanishes at infinity: 
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Equation ^ conserves H = f + U{X) + 2^(i^ + 
LiP'A^) + F{X)A, where the first two terms describe en- 
ergy in the propagating wave, the second two, energy in 
A and the final one, the coupling energy. 

The dynamics of X{t) in (|6|, neglecting A, conserves 
an energy E — ^X^ + U{X), and the trajectories he on 
level sets of E, figure |2] with E > along unbounded 
{E < along bounded) trajectories and separatrix or- 
bits along which E = 0. In the first instance the phase 




L(u, X, t) 



£{u{x, t)) dx dt. 



(5) 



ODEs are derived by assuming the solution de- 
pends on a few time-dependent parameters u(x,t) w 
Mansatz (-^1 (i) , • ■ • ,Xn(t)), inserting this ansatz into inte- 
gral ([5| and integrating out the x-dependence to obtain 
a finite-dimensional Lagrangian whose Euler-Lagrange 
equations describe the evolution of the parameters X(t). 



FIG. 2: The uncoupled phase portrait due to the potential 
U{X) in (a) equations ([7| and (b) equations ([8|. 



plane has two heteroclinic orbits connecting degenerate 
(saddle-type) fixed points at {X,X) — (±cxi,0), while 
the second has one homoclinic to (-Koo, 0). 

When A{t) is allowed to vary, the level sets of E 
cease to be invariant. Define capture (escape) to be a 
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trajectory that crosses a separatrix from the region of 
unbounded trajectories to that of bounded trajectories 
(bounded to unbounded) . We construct approximate so- 
lutions via matched asymptotic approximations where 
'outer solutions' consist of expansions near the degen- 
erate saddle points, which are connected via 'inner so- 
lutions,' i. e. separatrix orbits. An energy change calcu- 
lated over each separatrix orbit is used to match together 
two consecutive outer approximations near infinity. 

Over a full trajectory from one saddle-approach to the 
next, the total change in E is the Melnikov integral [5] 
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^ A^F{X{t))dt 
dt 
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F{X{t))—dt 



(9) 



assuming F{X{t)) ^ as \t\ oo. Equation ^ may 
be solved for A by variation of parameters and used to 
simplify ([9|. Under the assumption that A{t) — > as 
t — > — cxD, this is 



AE= 



F{Xs{t))e"^'dt 



(10) 



the Fourier transform of F{Xs{t)) evaluated at oj, the 
resonant frequency of A(t). Here X{t) has been approx- 
imated by Xgit), the solution along the separatrix. The 
critical velocity, which solves mv'^/2 = AE, is 



F{Xs{t))e 



(11) 



For system ([T]), Xg = sinh ^(t — ti) and Vc = tt-^i 
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(t-tiY 



and Vc — 



While for system ([8| , = log , , 

V2c7re^'^^'^^. These are shown in figure [s^ and b (with 
c = e^, TO = 1 and lj — e^^/^ in Q). 

The time ti is the 'symmetry time' of the first inner 
solution, at which Xs = in model ^ and -^Xs — 
in model If i? > after the first interaction, the 
wave moves off toward infinity. If i5 < the wave turns 
around. The oscillator A{t) has become excited and the 
variation of parameters formula shows that, for large t, 

A{t) -smuj{t~ti) [ F{Xs{t))cosuj{t -ti)dT 

^ J-oo 
C f°° 

H cosLu{t — ti) / F{Xs{t)) sin uj{t — ti)dT. 

^ J-oo 

(12) 

In model ([t]), F{Xs{t)) is odd about ti, so the first 
integral vanishes identically and A(t) oc cos ui(t — ti). 
In model ([s]), F{Xs{t)) is even about ti, so A{t) oc 
sin a;(t — ii)!^ The solution now alternates between two 
behaviors — negative energy 'outer solutions' dominated 
by the degenerate fixed point at ±oo and near-separatrix 
solutions with center time ti, until at some step n > 2, 
En > and the pulse escapes. 



Equation ( 10 1 gives AE along the first near-separatrix 
solution. At each subsequent interaction, at time ti, a 
similar calculation is performed [5J [3], with new terms 
that arise because A{t) no longer approaches zero as in 
backward time along the near-separatrix solution, but in- 
stead is asymptotically given by a sum of terms like ( 12 1, 



one for each previous collision. The energy level E after 
the nth interaction depends, thus, not only on the ini- 
tial energy, but on the sequence of times ti through t„. 
The time difference tj — tj-i, in turn, is a function of 
the energy level Ej_i, since the period of this nonlinear 
oscillator depends on its energy, as we show below. 

This time-change can be calculated by the matching 
conditions between the near-separatrix solution centered 
at tj_i and the near-saddle-expansion (outer solution) 
immediately following, and then the near saddle expan- 
sion to the next near-separatrix solution. 

Under the assumption that as X ^ oo, U{X) ~ 
—a'^e~^^^/2, we examine the large-X behavior of the 
first near-separatrix solution. A divergent integral for 
t — ti is regularized as: 
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A calculation for the ensuing large- X-saddle (outer) ap- 
proximate solution with energy Ei < Q, assumed to reach 
its maximum X = X* at t = t* , yields 



I-2E- 



cos 



-Pit* -t)] = 



V^2E: 



(14) 



Matching ( 14 1 with ( 13 1 yields, via a consistency con- 



dition, t* — ti. The calculation for t2 — t* is identical 
and eliminating t* gives t2 — ti . This calculation and its 
conclusion are unchanged for each time interval (<j-i, tj) 
and energy Ej-i, yielding 



tj — tj—i — 27^ 



-2E,^ip- 



(15) 



As t- ti 



oo along the jth near-separatrix solution. 



■J 2 /e 
A{t) ^Y.ak—ti-iguj{t~tk), 



(16) 



k=l 



where tifctrigx = (— 1)*^+"'^ cos a; in model ([t]) and 
CTfc trig X — ~ sin x in model . This implies 



-{l + 2j2'^j~kCOSLu{tj -tk)), (17) 



k=l 



with initial energy Eq = 
applied alternately, constitute the separatrix map. 



Equations (15 1 and (17 1 
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Summing equation (17), we find 



and find impressive agreement for e = 0.25. The reso- 



En — En — 
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(18) 



fe=l 3 = 1 



If En > 0, the waves escape to infinity with escape ve- 
locity Uout = '^\/'^En/m and t„+i undefined in (15). A 
solution is an n-bounce resonance if E^ — Eq, i.e. if 



n n 



= 2^ aj-k cosLoitk - tj) 

3 = 1 k=l 



71 n 



.j = l k=l 



i. e. (3„ = where Q„ = Yl"=i CTje"^*^ . 

For an exact n-bounce resonant solution, En-j = Ej. 
Letting 9j = Lu(t j^i — tj), this implies = 9j which 
yields the following conditions: 



If n = 2m 



Ifn = 2m + 1 



^(-1)^+1 sin (^0, 



m 

E 



■m 

(-1)^+1 cos (5^0.) 



= 0. 



For model ([t]), we find for 2- and 3-bounce resonances: 

„(2) _ 



"'^^ ^ r ' - 4(;rr^ 



(19) 



while for n > 4 these formulas must be solved numeri- 



cally, in conjunction with (15). 

In models (7| and ([s]), TZ — o(l) in equation (15), 
whereas TZ, calculated numerically, is 0(1) for the model 
studied in [B^. Thus (19) must be modified to include 
nonzero TZ. 

Most initial conditions do not lead to exact n-bounce 
resonances and we approximate the dynamics by iterat- 
ing the map defined by (15) and (17 1. The equivalent to 
figure [iji may be created by varying Ei^ = w?;fjj/2, and 



iterating until En = Eo 



t/2 > 0. In addition. 



the graph is color-coded by n, the number of bounces pre- 
ceding escape. In figures [Sj: and|3}i we compare solutions 
of this map to solutions of the sine-Gordon model ([t]) 



nant initial velocities given by equation (19) are marked 
in figure |3|:. In [B] we show, for the map describing ([T]), 
clusters of (n -I- l)-bounce windows accumulate at the 
edges of each n-bounce window, repeated at diminishing 
scales, with an intricacy that would be difficult to achieve 
from ODE initial- value simulations. For a wider view of 
fractal behavior, figure [S}? shows, using false-color, n, the 
number of bounces or interactions before the kink es- 
capes, as a function of both e and the initial velocity. 
The black curve gives Vc{e), and the fractal structure to 
its left shows how the windows appear and disappear as 
e is varied. 

In summary, we have explained the intricate chaotic 
dynamics arising in the interactions of solitary waves by 
the dynamics of simple iterated separatrix maps. Our 
method applies to many such systems studied over the 
past 25 years. We have shown that the 2-bounce reso- 
nance is the simplest manifestation of this chaotic scat- 
tering process. 
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FIG. 3: (Color Online) (a): The critical velocity Vc for 
model ([7|. (b): As (a), for model Q. (c): vin vs. Vout 
for the discrete map with the number of bounces coded by 
color; from ( |19[ ) vi^^ are marked by • and by T. (d): As 
(c), for numerical integration of ODE |[7|. (e): Number of 
bounces as a function of botli e and the initial velocity. 
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